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Abstract 



Mixture models are a fundamental tool in applied statistics and machine learning for treating 

>— I . data taken from multiple subpopulations. The current practice for estimating the parameters 

c/3 ' of such models relies on local search heuristics (e.g., the EM algorithm) which are prone to 

failure, and existing consistent methods are unfavorable due to their high computational and 

sample complexity which typically scale exponentially with the number of mixture components. 

This work develops an efficient method of moments approach to parameter estimation for a 

broad class of high-dimensional mixture models with many components, including multi-view 

mixtures of Gaussians (such as mixtures of axis-aligned Gaussians) and hidden Markov models. 

The new method leads to rigorous unsupervised learning results for mixture models that were 

not achieved by previous works; and, because of its simplicity, it offers a viable alternative to 

~^ ' EM for practical deployment. 
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1 Introduction 

Mixture models are a fundamental tool in applied statistics and machine learning for treating data 
taken from multiple subpopulations (Titterington et al., 1985). In a mixture model, the data are 
generated from a number of possible sources, and it is of interest to identify the nature of the 
individual sources. As such, estimating the unknown parameters of the mixture model from sam- 
pled data — especially the parameters of the underlying constituent distributions — is an important 
statistical task. For most mixture models, including the widely used mixtures of Gaussians and 
hidden Markov models (HMMs), the current practice relies on the Expectation-Maximization (EM) 
algorithm, a local search heuristic for maximum likelihood estimation. However, EM has a number 
of well-documented drawbacks regularly faced by practitioners, including slow convergence and 
suboptimal local optima (Redner and Walker, 1984). 

An alternative to maximum likelihood and EM, especially in the context of mixture models, 
is the method of moments approach. The method of moments dates back to the origins of mix- 
ture models with Pearson's solution for identifying the parameters of a mixture of two univariate 
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Gaussians (Pearson, 1894). In this approach, model parameters are chosen to specify a distribution 
whose p-th. order moments, for several values of p, are equal to the corresponding empirical mo- 
ments observed in the data. Since Pearson's work, the method of moments has been studied and 
adapted for a variety of problems; their intuitive appeal is also complemented with a guarantee of 
statistical consistency under mild conditions. Unfortunately, the method often runs into trouble 
with large mixtures of high-dimensional distributions. This is because the equations determining 
the parameters are typically based on moments of order equal to the number of model parameters, 
and high-order moments are exceedingly difficult to estimate accurately due to their large variance. 
This work develops a computationally efficient method of moments based on only low-order 
moments that can be used to estimate the parameters of a broad class of high-dimensional mixture 
models with many components. The resulting estimators can be implemented with standard nu- 
merical linear algebra routines (singular value and eigenvalue decompositions), and the estimates 
have low variance because they only involve low-order moments. The class of models covered by 
the method includes certain multivariate Gaussian mixture models and HMMs, as well as mix- 
ture models with no explicit likelihood equations. The method exploits the availability of multiple 
indirect "views" of a model's underlying latent variable that determines the source distribution, 
although the notion of a "view" is rather general. For instance, in an HMM, the past, present, 
and future observations can be thought of as different noisy views of the present hidden state; in 
a mixture of product distributions (such as axis-aligned Gaussians) , the coordinates in the output 
space can be partitioned (say, randomly) into multiple non-redundant "views". The new method 
of moments leads to unsupervised learning guarantees for mixture models under mild rank condi- 
tions that were not achieved by previous works; in particular, the sample complexity of accurate 
parameter estimation is shown to be polynomial in the number of mixture components and other 
relevant quantities. Finally, due to its simplicity, the new method (or variants thereof) also offers 
a viable alternative to EM and maximum likelihood for practical deployment. 

1.1 Related work 

Gaussian mixture models. The statistical literature on mixture models is vast (a more thorough 
treatment can be found in the texts of Titterington et al. (1985) and Lindsay (1995)), and many 
advances have been made in computer science and machine learning over the past decade or so, in 
part due to their importance in modern applications. The use of mixture models for clustering data 
comprises a large part of this work, beginning with the work of Dasgupta (1999) on learning mix- 
tures of k well-separated (i-dimensional Gaussians. This and subsequent work (Arora and Kannan, 
2001; Dasgupta and Schulman, 2007; Vempala and Wang, 2002; Kannan et al., 2005; Achlioptas 
and McSherry, 2005; Chaudhuri and Rao, 2008; Brubaker and Vempala, 2008; Chaudhuri et al., 
2009) have focused on efficient algorithms that provably recover the parameters of the constituent 
Gaussians from data generated by such a mixture distribution, provided that the distance between 
each pair of means is sufficiently large (roughly either d c or k c times the standard deviation of the 
Gaussians, for some c > 0). Such separation conditions are natural to expect in many clustering 
applications, and a number of spectral projection techniques have been shown to enhance the sep- 
aration (Vempala and Wang, 2002; Kannan et al., 2005; Brubaker and Vempala, 2008; Chaudhuri 
et al., 2009). More recently, techniques have been developed for learning mixtures of Gaussians 
without any separation condition (Kalai et al., 2010; Belkin and Sinha, 2010; Moitra and Valiant, 
2010), although the computational and sample complexities of these methods grow exponentially 
with the number of mixture components k. This dependence has also been shown to be inevitable 



without further assumptions (Moitra and Valiant, 2010). 

Method of moments. The latter works of Belkin and Sinha (2010), Kalai et al. (2010), and Moitra 
and Valiant (2010) (as well as the algorithms of Feldman et al. (2005, 2006) for a related but different 
learning objective) can be thought of as modern implementations of the method of moments, and 
their exponential dependence on k is not surprising given the literature on other moment methods 
for mixture models. In particular, a number of moment methods for both discrete and continuous 
mixture models have been developed using techniques such as the Vandermonde decompositions 
of Hankel matrices (Lindsay, 1989; Lindsay and Basak, 1993; Boley et al., 1997; Gravin et al., 
2012). In these methods, following the spirit of Pearson's original solution, the model parameters 
are derived from the roots of polynomials whose coefficients are based on moments up to the 
il(/c)-th order. The accurate estimation of such moments generally has computational and sample 
complexity exponential in k. 

Spectral approach to parameter estimation with low-order moments. The present work 
is based on a notable exception to the above situation, namely Chang's spectral decomposition 
technique for discrete Markov models of evolution (Chang, 1996) (see also Mossel and Roch (2006) 
and Hsu et al. (2009) for adaptations to other discrete mixture models such as discrete HMMs). 
This spectral technique depends only on moments up to the third-order; consequently, the resulting 
algorithms have computational and sample complexity that scales only polynomially in the number 
of mixture components k. The success of the technique depends on a certain rank condition of 
the transition matrices; but this condition is much milder than separation conditions of clustering 
works, and it remains sufficient even when the dimension of the observation space is very large (Hsu 
et al., 2009). In this work, we extend Chang's spectral technique to develop a general method of 
moments approach to parameter estimation, which is applicable to a large class of mixture models 
and HMMs with both discrete and continuous component distributions in high-dimensional spaces. 
Like the moment methods of Moitra and Valiant (2010) and Belkin and Sinha (2010), our algorithm 
does not require a separation condition; but unlike those previous methods, the algorithm has 
computational and sample complexity polynomial in k. 

Some previous spectral approaches for related learning problems only use second-order moments, 
but these approaches can only estimate a subspace containing the parameter vectors and not the 
parameters themselves (McSherry, 2001). Indeed, it is known that the parameters of even very sim- 
ple discrete mixture models are not generally identifiable from only second-order moments (Chang, 
1996) 1 . We note that moments beyond the second-order (specifically, fourth-order moments) have 
been exploited in the methods of Frieze et al. (1996) and Nguyen and Regev (2009) for the prob- 
lem of learning a parallelepiped from random samples, and that these methods are very related to 
techniques used for independent component analysis (Hyvarinen and Oja, 2000). Adapting these 
techniques for other parameter estimation problems is an enticing possibility. 

Multi-view learning. The spectral technique we employ depends on the availability of multiple 
views, and such a multi-view assumption has been exploited in previous works on learning mixtures 
of well-separated distributions (Chaudhuri and Rao, 2008; Chaudhuri et al., 2009). In these previous 
works, a projection based on a canonical correlation analysis (Hotelling, 1935) between two views is 
used to reinforce the separation between the mixture components, and to cancel out noise orthogonal 
to the separation directions. The present work, which uses similar correlation-based projections, 
shows that the availability of a third view of the data can remove the separation condition entirely. 



See Appendix G for an example of Chang (1996) demonstrating the non-identifiability of parameters from only 
second-order moments in a simple class of Markov models. 



The multi-view assumption substantially generalizes the case where the component distributions 
are product distributions (such as axis-aligned Gaussians), which has been previously studied in the 
literature (Dasgupta, 1999; Vempala and Wang, 2002; Chaudhuri and Rao, 2008; Feldman et al., 
2005, 2006); the combination of this and a non-degeneracy assumption is what allows us to avoid 
the sample complexity lower bound of Moitra and Valiant (2010) for Gaussian mixture models. The 
multi-view assumption also naturally arises in many applications, such as in multimedia data with 
(say) text, audio, and video components (Blaschko and Lampert, 2008; Chaudhuri et al., 2009); as 
well as in linguistic data, where the different words in a sentence or paragraph are considered noisy 
predictors of the underlying semantics (Gale et al., 1992). In the vein of this latter example, we 
consider estimation in a simple bag-of-words document topic model as a warm-up to our general 
method; even this simpler model illustrates the power of pair-wise and triple- wise (i.e., bigram and 
trigram) statistics that were not exploited by previous works on multi-view learning. 

1.2 Outline 

Section 2 first develops the method of moments in the context of a simple discrete mixture model 
motivated by document topic modeling; an explicit algorithm and convergence analysis are also 
provided. The general setting is considered in Section 3, where the main algorithm and its accom- 
panying correctness and efficiency guarantee are presented. Applications to learning multi-view 
mixtures of Gaussians and HMMs are discussed in Section 4. All proofs are given in the appendix. 

1.3 Notations 

The standard inner product between vectors u and v is denoted by (u, v) = u T v. We denote the 
p-norm of a vector vby \\v\\ p . For a matrix A G M mxn , we let \\A\\2 denote its spectral norm ||.A|| 2 := 
sup-/(j||.Au||2/||w||2, ||-A||f denote its Frobenius norm, (Ji(A) denote the i-th. largest singular value, 
and k(A) := cri(A)/a m i n r m ^ n \(A) denote its condition number. Let A n_1 := {(pi,P2, ■ ■ ■ ,Pn) G R n : 
Pi > Vi, Y17=iPi = 1} denote the probability simplex in M. n , and let S n ~ l := {u G W 1 : ||u||2 = 1} 
denote the unit sphere in M. n . Let e*j G K denote the i-th coordinate vector whose i-th entry is 1 
and the rest are zero. Finally, for a positive integer n, let [n] := {1, 2, . . . , n}. 

2 Warm-up: bag-of-words document topic modeling 

We first describe our method of moments in the simpler context of bag-of-words models for docu- 
ments. 

2.1 Setting 

Suppose a document corpus can be partitioned by topic, with each document being assigned a single 
topic. Further, suppose the words in a document are drawn independently from a multinomial 
distribution corresponding to the document's topic. Let k be the number of distinct topics in the 
corpus, d be the number of distinct words in the vocabulary, and £ > 3 be the number of words in 
each document (so the documents may be quite short). 

The generative process for a document is given as follows: 



1. The document's topic is drawn according to the multinomial distribution specified by the 
probability vector w = (wi,W2, ■ ■ ■ ,Wk) G A . This is modeled as a discrete random 
variable h such that 

Pr[/i = j] = Wj , j G [k]. 

2. Given the topic h, the document's £ words are drawn independently according to the multi- 
nomial distribution specified by the probability vector p,^ G A d_1 . The random vectors 
x\,X2, ■ ■ ■ ,xe £M. d represent the £ words by setting 

x v = e i 44> the u-th word in the document is i, i & [d] 

(the reason for this encoding of words will become clear in the next section) . Therefore, for 
each word v G [£] in the document, 

Pr[x v = &i\h = j] = {ei,flj) = Mi j, i G [d],j G [k], 

where M G M. dxk is the matrix of conditional probabilities M := \jl1\jl2\ ■ ■ ■ \fik\- 

This probabilistic model has the conditional independence structure depicted in Figure 2(a) as a 
directed graphical model. 

We assume the following condition on w and M. 

Condition 2.1 (Non-degeneracy: document topic model). Wj >0 for all j€ [k], and M has rank k. 

This condition requires that each topic has non-zero probability, and also prevents any topic's 
word distribution from being a mixture of the other topics' word distributions. 

2.2 Pair-wise and triple-wise probabilities 

Define Pairs G M. dxd to be the matrix of pair-wise probabilities whose (i,j)-th entry is 

Pairsjj := Pr[a?i = e*i,f 2 = e}], i,j G [d]. 

Also define Triples G ]R" X " X " to be the third-order tensor of triple- wise probabilities whose (i,j, k)- 
th entry is 

TripleSjj K := Pr[xi = e*, x 2 = ej,x 3 = e K ], i, j, k G [d\. 

The identification of words with coordinate vectors allows Pairs and Triples to be viewed as expec- 
tations of tensor products of the random vectors xi, X2, and X3: 

Pairs = E[xi C3 X2] and Triples = E[xi X2 <8> S3]. (1) 

We may also view Triples as a linear operator Triples: R — > R given by 

Triples(^) := E[(x\ <g> X2)(ff,x 3 )]. 

In other words, the (i,j)-th entry of Triples(??) for ff = (r]i,r]2, ■ ■ ■ ,%) is 

d d 

Triples(^)ij = ^ ^Triples^ ? - x = ^r/ x Triples(e x )i,j- 

x=l x=l 

The following lemma shows that Pairs and Triples (ff) can be viewed as certain matrix products 
involving the model parameters M and w. 



Lemma 2.1. Pairs = Mdiag(u>)M T and Triples^) = Mdiag(M T r/) di&g(w)M T for all ff 6 
Proof. Since a?i, X2, and x% are conditionally independent given h, 

k 

Pairsjj = Pr[xi = e*j, x 2 = ej] = y, Pr[^i = ej, x 2 = &j\h = t] ■ Pi[h = t] 

t=i 

k k 

= ^2 Pr ^i = e,j\h = t\- Pr[x 2 = Bj\h = t] ■ Pr[/i = t] = ^ M M • M ijt • w t 
t=i t=\ 

so Pairs = M diag(uf)M T . Moreover, writing ff = (771, 772, • • • , %), 

d 



Triples(f/)jj = ^^11 Pr[f 1 = e h x 2 = ej,x 3 = e x ] 

x=l 

d k k 

= Yl Yl Vx ■ M i,t ■ M jyt ■ M Xyt ■w t = J2 M M • M jyt ■ w t ■ (M T ff) t 



x=l t=\ t=\ 

so Triples^) = M diag(M T ff) diag(tt;)M T . D 

2.3 Observable operators and their spectral properties 

The pair-wise and triple-wise probabilities can be related in a way that essentially reveals the con- 
ditional probability matrix M. This is achieved through a matrix called an "observable operator". 
Similar observable operators were previously used to characterize multiplicity automata (Schiitzenberger, 
1961; Jaeger, 2000) and, more recently, for learning discrete HMMs (via an operator parameteriza- 
tion) (Hsu et al., 2009). 

Lemma 2.2. Assume Condition 2.1. Let U G R dxfc and V G M. dxk be matrices such that both 
U T M and V T M are invertible. Then [7 T PairsV r is invertible, and for all ff £ R rf , the "observable 
operator" B{ff) £ R , given by 

B(ff) := (?7 T Triples(??)y)(c/ T Pairsy)- 1 , 

satisfies 

B(ff) = (U T M)diag(M T fj){U T My 1 . 

Proof. Since diag(u>) >- by Condition 2.1 and [7 T PairsV = (U T M) diag(w)M T V by Lemma 2.1, it 
follows that U T PairsV is invertible by the assumptions on U and V. Moreover, also by Lemma 2.1, 

B(fj) = (^ T Triples(7?)y) (c/ T PairsV r )" 1 

= (?7 T Mdiag(M T ^)diag(u;)M T T/) (c^PairsV)- 1 

= {U T M)diag{M T ff)(U T M)- 1 (c/ T Mdiag(u;)M T V) (l^PairsV)- 1 

= (U T M)diag{M T ff)(U T M)- 1 . □ 

The matrix B{rf) is called "observable" because it is only a function of the observable variables' 
joint probabilities {e.g., Pr[a?i = ei,x 2 = ej]). In the case ff = e x for some x £ [d], the matrix 

6 



Algorithm A 

1. Obtain empirical frequencies of word pairs and triples from a given sample of documents, and 
form the tables Pairs <G R dxd and Triples <G M. dxdxd corresponding to the population quantities 
Pairs and Triples. 

2. Let U £ R dxfc and V £ M. dxk be, respectively, matrices of orthonormal left and right singular 
vectors of Pairs corresponding to its top k singular values. 

3. Pick fj £ M. d (see remark in the main text), and compute the right eigenvectors £i, £2, • • • , £fc (of 
unit Euclidean norm) of 

B(fj) := (U T TT\p\cs(rj)V)(U T P^sV ) _1 . 
(Fail if not possible.) 

4. Let fj,j := u£ 3 /(T, tlij) for all 3 £ [k]. 

5. Return M :— [/&1I/&2I ■ ■ ■ |Afc]- 



Figure 1: Topic- word distribution estimator (Algorithm A). 

B{e x ) is similar (in the linear algebraic sense) to the diagonal matrix diag(M T e x ); the collection 
of matrices {diag(M T e* r ) : x £ [d]} (together with w) can be used to compute joint probabilities 
under the model (see, e.g., Hsu et al. (2009)). Note that the columns of U T M are eigenvectors of 
B(e x ), with the j'-th column having an associated eigenvalue equal to Pr[x„ = x\h = j]. If the word 
x has distinct probabilities under every topic, then B(e x ) has exactly k distinct eigenvalues, each 
having geometric multiplicity one and corresponding to a column of U T M. 

2.4 Topic-word distribution estimator and convergence guarantee 

The spectral properties of the observable operators B(fj) implied by Lemma 2.2 suggest the estima- 
tion procedure (Algorithm A) in Figure 1. The procedure is essentially a plug-in approach based on 
the equations relating the second- and third-order moments in Lemma 2.2. We focus on estimating 
M; estimating the mixing weights w is easily handled as a secondary step (see Appendix B.5 for 
the estimator in the context of the general model in Section 3.1). 

On the choice of fj. As discussed in the previous section, a suitable choice for fj can be based 
on prior knowledge about the topic-word distributions, such as fj = e x for some x £ [d] that has 
different conditional probabilities under each topic. In the absence of such information, one may 
select fj randomly from the subspace range(£7). Specifically, take fj := U0 where 9 £ M. k is a random 
unit vector distributed uniformly over S . 

The following theorem establishes the convergence rate of Algorithm A. 

Theorem 2.1. There exists a constant C > such that the following holds. Pick any 5 £ (0, 1). 
Assume the document topic model from Section 2.1 satisfies Condition 2.1. Further, assume that 
in Algorithm A, Pairs and Triples are, respectively, the empirical averages of N independent copies 
of x\ ® X2 and x\ <£> X2 ® X3 ; and that fj = U9 where 9 £ R is an independent random unit vector 
distributed uniformly over S k ~ 1 . If 

n>c tMn(1/a) 



CTfc(M) 6 -(Tfc (Pairs) 4 -<5 2 ' 






> x 2 

(aj (b) 

Figure 2: (a) The multi-view mixture model, (b) A hidden Markov model. 

then with probability at least 1 — 5, the parameters returned by Algorithm A have the following 

guarantee: there exists a permutation r on [k] and scalars ci,C2, . . . , c& G M such that, for each 

3 6 [fe], 

II - - II <r II- II fc 5 / jgVg 

||c,-Mi -Mr(,-)»2 < C- llMrwIla • ^4.^(^)2., • V ^V^- 

The proof of Theorem 2.1, as well as some illustrative empirical results on using Algorithm A, 
are presented in Appendix A. A few remarks about the theorem are in order. 

On boosting the confidence. Although the convergence depends polynomially on 1/(5, where 5 
is the failure probability, it is possible to boost the confidence by repeating Step 3 of Algorithm 
A with different random fj until the eigenvalues of B(ff) are sufficiently separated (as judged by 
confidence intervals). 

On the scaling factors Cj. With a larger sample complexity that depends on d, an error bound 
can be established for \\p,j — fJL r (j\\\i directly (without the unknown scaling factors Cj). We also 
remark that the scaling factors can be estimated from the eigenvalues of B(ff), but we do not 
pursue this approach as it is subsumed by Algorithm B anyway. 

3 A method of moments for multi-view mixture models 

We now consider a much broader class of mixture models and present a general method of moments 
in this context. 

3.1 General setting 

Consider the following multi-view mixture model; k denotes the number of mixture components, 
and £ denotes the number of views. We assume t > 3 throughout. Let w = {w\,W2, ■ ■ ■ , Wk) £ A fc_1 
be a vector of mixing weights, and let h be a (hidden) discrete random variable with Pr[/i = j] = Wj 
for all j G [A;]. Let a?i, X2, ■ ■ ■ , xe £ K d be f random vectors that are conditionally independent given 
h; the directed graphical model is depicted in Figure 2(a). 
Define the conditional mean vectors as 

j2 v>j := E[x v \h = j], v G [£],j G [k], 

and let M v G M. dxk be the matrix whose j'-th column is j3 v ,j- Note that we do not specify anything 
else about the (conditional) distribution of x v — it may be continuous, discrete, or even a hybrid 
depending on h. 

We assume the following conditions on w and the M v . 



Condition 3.1 (Non-degeneracy: general setting). Wj > for all j G [k], and M v has rank k for 
alive [£]. 

We remark that it is easy to generalize to the case where views have different dimensionality 
(e.g., x v G M. dv for possibly different dimensions d v ). For notational simplicity, we stick to the same 
dimension for each view. Moreover, Condition 3.1 can be relaxed in some cases; we discuss one 
such case in Section 4.1 in the context of Gaussian mixture models. 

Because the conditional distribution of x v is not specified beyond its conditional means, it is 
not possible to develop a maximum likelihood approach to parameter estimation. Instead, as in the 
document topic model, we develop a method of moments based on solving polynomial equations 
arising from eigenvalue problems. 

3.2 Observable moments and operators 

We focus on the moments concerning {x\,X2,x^}, but the same properties hold for other triples of 
the random vectors {x a ,Xb,x c } C {x v : v G [£]} as well. 

As in (1), we define the matrix Pi, 2 G M rfxd of second-order moments, and the tensor -Pi, 2, 3 G 
M of third-order moments, by 

Pi, 2 := E[a?i <g> x 2 ] and Pi,2,3 : = E [^i ®x 2 ® x 3 ]. 

Again, Pi,2,3 is regarded as the linear operator Pi.,2,3: ff ^ E[(a?i <8) X2){ff, X3)]. 

Lemma 3.1 and Lemma 3.2 are straightforward generalizations of Lemma 2.1 and Lemma 2.2. 

Lemma 3.1. P l ^ = M 1 diag(w)M 2 T and P\,2,?>{ff) = M 1 diag(M 3 T ?7) diag(w;)M 2 T for all ff G R d . 

Lemma 3.2. Assume Condition 3.1. For v G {1,2,3}, let U v G M. dxk be a matrix such that 
Uj M v is invertible. Then U1P12U2 is invertible, and for all ff G M , the "observable operator" 
Bi,2,s(v) 6 RkXk > 9iven by B 1>2 ,M := (U^ P^Mv^)^ Pi^)' 1 , satisfies 

Bi,2,M = (f/ 1 T M 1 )diag(M 3 T r/)(C/ 1 T M 1 )- 1 . 
In particular, the k roots of the polynomial A i-> det(i?i ) 2,3(f/) — A/) are {{ff,p3j) '■ j G [fc]}. 

Recall that Algorithm A relates the eigenvectors of B(ff) to the matrix of conditional means M. 
However, eigenvectors are only defined up to a scaling of each vector; without prior knowledge of 
the correct scaling, the eigenvectors are not sufficient to recover the parameters M. Nevertheless, 
the eigenvalues also carry information about the parameters, as shown in Lemma 3.2, and it is 
possible to reconstruct the parameters from different the observation operators applied to different 
vectors ff. This idea is captured in the following lemma. 

Lemma 3.3. Consider the setting and definitions from Lemma 3.2. Let G K fcxfc be an invertible 
matrix, and let QJ G M fc be its i-th row. Moreover, for all i G [A;], let A^i, Aj^, . . . , \i,k denote the 
k eigenvalues of -61,2,3(^3^1) *n the order specified by the matrix of right eigenvectors U^Mi. Let 
L G R be the matrix whose (i,j)-th entry is Ajj. Then 

QUJM-i = L. 



Observe that the unknown parameters M3 are expressed as the solution to a linear system in 
the above equation, where the elements of the right-hand side L are the roots of fe-th degree poly- 
nomials derived from the second- and third-order observable moments (namely, the characteristic 
polynomials of the -61,2,3(^3^1)) Vi £ [k]). This template is also found in other moment methods 
based on decompositions of a Hankel matrix. A crucial distinction, however, is that the k-ih degree 
polynomials in Lemma 3.3 only involve low-order moments, whereas standard methods may involve 
up to fi(/s)-th order moments which are difficult to estimate (Lindsay, 1989; Lindsay and Basak, 
1993; Gravin et al., 2012). 

3.3 Main result: general estimation procedure and sample complexity bound 

The lemmas in the previous section suggest the estimation procedure (Algorithm B) presented in 
Figure 3. 



Algorithm B 

1. Compute empirical averages from N independent copies of X\®X2 to form Pi, 2 £ K dxd . Similarly 
do the same for x\ £3 to form P13 £ M. kxk , and for x 1 $2 x$ to form A, 2, 3 G R dxdxd . 

2. Let U\ £ W lxk and U2 £ R dxfc be, respectively, matrices of orthonormal left and right singular 
vectors of Pi, 2 corresponding to its top k singular values. Let C/3 £ R dxk be the matrix of 
orthonormal right singular vectors of Pi, 3 corresponding to its top k singular values. 

3. Pick an invertible matrix 6 £ R fcxfc , with its z-th row denoted as 8J £ M. k . In the absence of any 
prior information about M3, a suitable choice for is a random rotation matrix. 

Form the matrix ^1,2,3(^3^1) := (#rA,2,3(#3&)#2)(#rA,2#2) _1 . 

Compute Pi £ ~R kxk (with unit Euclidean norm columns) that diagonalizes -81,2,3(^3^1), i.e., 

Ri 1 Bi : 2, z{UzQi)Ri = diag(Ai,i, Ai,2, . . . , Ai,fc). (Fail if not possible.) 

4. For each i £ {2, . . . , k}, obtain the diagonal entries Aj,i, A»,2, • • • , K,k of P 1 _1 -Bi,2,3(^3^i)Pi, and 
form the matrix L £ R kxk whose (i,j)-ih entry is Xij. 

5. Return M 3 := U^Q^L. 



Figure 3: General method of moments estimator (Algorithm B). 

As stated, the Algorithm B yields an estimator for M3, but the method can easily be applied 
to estimate M v for all other views v. One caveat is that the estimators may not yield the same 
ordering of the columns, due to the unspecified order of the eigenvectors obtained in the third step 
of the method, and therefore some care is needed to obtain a consistent ordering. We outline one 
solution in Appendix B.4. 

The sample complexity of Algorithm B depends on the specific concentration properties of 
xi,X2,x*3- We abstract away this dependence in the following condition. 

Condition 3.2. There exist positive scalars Nq, C12, Ci,3, Ci,2,3> and a function f(N, 5) (decreas- 
ing in A^ and S) such that for any N > Nq and 6 £ (0, 1), 



1. Pr 



P a ,b-Pa,bh<C a ,b-f(N,S) > 1 - S for {a, b} € {{1, 2}, {1, 3}}, 



2. W e R d , Pr 



1^1,2,3 (0) " ^1,2,3^)112 < Ci, 2 ,3 ■ N| 2 ■ f(N, 5) 
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> 1 - s. 



Moreover (for technical convenience), Pi 3 is independent of Pi, 2, 3 (which may be achieved, say, by 
splitting a sample of size 27V). 

For the discrete models such as the document topic model of Section 2.1 and discrete HMMs (Mos- 
sel and Roch, 2006; Hsu et al., 2009), Condition 3.2 holds with N = C 1>2 = Ci,3 = Ci, 2 ,3 = 1, and 
f(N,5) = (l + y / ln(l/<5))/v / 7V- Using standard techniques (e.g., Chaudhuri et al. (2009); Vershynin 
(2012)), the condition can also be shown to hold for mixtures of various continuous distributions 
such as multivariate Gaussians. 

Now we are ready to present the main theorem of this section (proved in Appendix B.6). 

Theorem 3.1. There exists a constant C > such that the following holds. Assume the three- 
view mixture model satisfies Condition 3.1 and Condition 3.2. Pick any e G (0, 1) and 5 G (0, Sq). 
Further, assume G G M. kxk is an independent random rotation matrix distributed uniformly over 
the Stiefel manifold {Q G M fcxfc : Q T Q = I}. If the number of samples N satisfies N > Nq and 

/( ' /k) - c ll2 ,3 • fc 5 • K{ Ml y Hk/6) ' e ' 



f(N,5) < C-min< 



min^j ||M 3 (ei - e})|| 2 • o- fc (Pi, 2 ) $ a k (P ] 



1,3; 



C\, 2 ■ \\Pi,2,3h ■ k^ ■ k(M 1 )^ hi(k/6y C 1;3 



where ||Pi, 2,31b := inax ^o 11-^1,2,3(^)112; then with probability at least 1 — 55, Algorithm B returns 

-^3 = [/i.3,1 1 A3,2| • • • |M3,ifc] with the following guarantee: there exists a permutation r on [k] such that 
for each j € [k] , 



\h,j ~ M3,r(j)l|2 < max p 3 ,j'||2 • e- 



4 Applications 



In addition to the document clustering model from Section 2, a number of natural latent variable 
models fit into this multi-view framework. We describe two such cases in this section: Gaussian 
mixture models and HMMs, both of which have been (at least partially) studied in the literature. 
In both cases, the estimation technique of Algorithm B leads to new learnability results that were 
not achieved by previous works. 

4.1 Multi-view Gaussian mixture models 

The standard Gaussian mixture model is parameterized by a mixing weight Wj, mean vector 
flj G 1H D , and covariance matrix Uj G ^L DxD for each mixture component j G [k]. The hidden 
discrete random variable h selects a component j with probability Pr[/i = j] = Wj] the conditional 
distribution of the observed random vector x given h is a multivariate Gaussian with mean /Z/j and 
covariance E^- 

The multi-view assumption for Gaussian mixture models asserts that for each component j, 
the covariance Uj has a block diagonal structure Uj = brkdiag(Zij, U 2 j, ■ ■ ■ , U^j) (a special case 
is an axis-aligned Gaussian). The various blocks correspond to the £ different views of the data 
xi, x 2 , . . . , X£ G M. d (for d = D/£), which are conditionally independent given h. The mean vector 
for each component j is similarly partitioned into the views as flj = (f2ij,i2 2 j, ■ ■ ■ ,fle,j)- Note that 
in the case of an axis-aligned Gaussian, each covariance matrix Uj is diagonal, and therefore the 
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original coordinates [D] can be partitioned into £ = 0(D/k) views (each of dimension d = Q(k)) in 
any way (say, randomly) provided that Condition 3.1 holds. 

Condition 3.1 requires that the conditional mean matrix M v = [(2 v ,i\j2 v ,2\ ■ ■ ■ \flv,k] f° r each view 
v have full column rank. This is similar to the non-degeneracy and spreading conditions used in 
previous studies of multi-view clustering (Chaudhuri and Rao, 2008; Chaudhuri et al., 2009). In 
these previous works, the multi-view and non-degeneracy assumptions are shown to reduce the 
minimum separation required for various efficient algorithms to learn the model parameters. In 
comparison, Algorithm B does not require a minimum separation condition at all. See Appendix D.3 
for details. 

While Algorithm B recovers just the means of the mixture components (see Appendix D.4 for 
details concerning Condition 3.2), we remark that a slight variation can be used to recover the 
covariances as well. Note that 

E[x v ®x v \h] = (M v e h ) (g> (M v e h ) + E Vjh = jl Vth ® jl Vih + E Vjh 

for all v G [£]. For a pair of vectors (f> £ M. d and ip G M. d , define the matrix Qi,2,3(<l>,4>) G M. dxd of 
fourth-order moments by Qi,2,3(</>, VO := E[(xi ^ X2){(p,X3)(ip,Xs}]. 

Proposition 4.1. Under the setting of Lemma 3.2, the matrix given by 

satisfies F lt2 ,3(4>,ip) = (C7' 1 r Mi)diag((0,//3 j t)(V',M3,t) + (<P,E 3tt ip) : t G [k])(U^ M^ 1 and hence is 
diagonalizable (in fact, by the same matrices as B\^,zijf))- 

Finally, we note that even if Condition 3.1 does not hold {e.g., if j2 v j = m G M rf (say) for all 
v € [£] , j G [k] so all of the Gaussians have the same mean) , one may still apply Algorithm B to 
the model (h,yi,ij2, ■ ■ ■ ,ye) where y v G m d + d ( d + 1 )/ 2 [ s the random vector that include both first- 
and second-order terms of x v , i.e., y v is the concatenation of x v and the upper triangular part of 
x v C5 x v . In this case, Condition 3.1 is replaced by a requirement that the matrices 



M' v := [ E[y v \h = 1] | E[y v \h = 2] | ■ ■ ■ | E[y v \h = k] ] G 



)(d+d(d+l)/2)xk 



of conditional means and covariances have full rank. This requirement can be met even if the means 
jl v j of the mixture components are all the same. 

4.2 Hidden Markov models 

A hidden Markov model is a latent variable model in which a hidden state sequence hi, h.2, ■ ■ ■ , hi 
forms a Markov chain hi — > h 2 — >■ ■••—>■ he over k possible states [k]; and given the state ht at 
time t G [k], the observation xt at time t (a random vector taking values in M. d ) is conditionally 
independent of all other observations and states. The directed graphical model is depicted in 
Figure 2(b). 

The vector n G A is the initial state distribution: 

Pr[hi = i] = iTi, i G [A;]. 
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For simplicity, we only consider time-homogeneous HMMs, although it is possible to generalize to 
the time- varying setting. The matrix T G R is a stochastic matrix describing the hidden state 
Markov chain: 

Pi[ht+i = i\h t = j] = T id , i,j G [k],t e[£- 1]. 

Finally, the columns of the matrix O = [o"i 1 02 1 • • • |ofc] £ R are the conditional means of the 
observation x t at time t given the corresponding hidden state h t : 

E[x t \h t = i] = Oei = di, i€[k],t€[£\. 

Note that both discrete and continuous observations are readily handled in this framework. For 
instance, the conditional distribution of Xt given ht = i (for i G [k]) could be a high-dimensional mul- 
tivariate Gaussian with mean 3i G R rf . Such models were not handled by previous methods (Chang, 
1996; Mossel and Roch, 2006; Hsu et al., 2009). 

The restriction of the HMM to three time steps, say t G {1,2,3}, is an instance of the three- view 
mixture model. 

Proposition 4.2. If the hidden variable h (from the three-view mixture model of Section 3.1) is 
identified with the second hidden state hi, then {x\,X2,x^} are conditionally independent given h, 
and the parameters of the resulting three-view mixture model on {h,x\,X2,x 3 ) are 

w := Ttt, Mi := Odiag(vf)T T diag(Tvf)- 1 , M 2 := O, M 3 := OT. 

From Proposition 4.2, it is easy to verify that -63,1, 2(77) = (UJOT) diag(0 T ff)(U 3 OT)^ 1 . There- 
fore, after recovering the observation conditional mean matrix O using Algorithm B, the Markov 
chain transition matrix can be recovered using the matrix of right eigenvectors R of .B3 1 2(77) and 
the equation (U^ 0)~ l R = T (up to scaling of the columns). 
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A Analysis of Algorithm A 

In this appendix, we give an analysis of Algorithm A (but defer most perturbation arguments to 
Appendix C), and also present some illustrative empirical results on text data using a modified 
implementation. 
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A.l Accuracy of moment estimates 



Lemma A.l. Fix 5 £ (0, 1). Let Pairs be the empirical average of N independent copies of x\®xi, 
and let Triples be the empirical average of N independent copies of {x\ ® X2){fj,x^). Then 



1. Pr 



2. Pr 



| Pairs — Pairs||F < 



1 + VM1/5) 



N 



> 1 — 5, and 



Mfj e R d , 1 1 Triples (ff) - Triples (rf)|| F < 



\fjui + VHU$J) 



N 



> 1-5. 



Proof. The first claim follows from applying Lemma F.l to the vectorizations of Pairs and Pairs 
(whereupon the Frobenius norm is the Euclidean norm of the vectorized matrices). For the second 
claim, we also apply Lemma F.l to Triples and Triples in the same way to obtain, with probability 
at least 1 — 5, 

E E B^W - Triples,, J 2 < (1 + V^W 

i=l j=l x=l 

Now condition on this event. For any fj = (j]i, 772, . . . , %) G M rf , 



^^(Triplesjj^ - Triples,^) 



x=l 



d d 

||Triples(r/) - Triples (r7)|| F = EE 

i=i i=i 

d d d 

^ E E 11 ^11 2 E^^^j^ ~ Tri P les i,i,x 

i=l j=l x=l 



< 



N 



where the first inequality follows by Cauchy-Schwarz. 



□ 



A. 2 Proof of Theorem 2.1 

Let E\ be the event in which 



I Pairs — Pairs llo < 



i + y/W/s) 



N 



and 



I J Triples (v) — Triples (v)\\ 2 < 



| 2 (i + vMV*)) 



N 



(3) 



for all v £ ~M. d . By Lemma A.l, a union bound, and the fact that ||A||2 < ||^4||f> we have Pr[£i] > 
1 — 25. Now condition on E\, and let E2 be the event in which 



7:=min|([/#,M(e;-e j ))| = min|(fl, [/ T M(e; - e}))| > 



V2a k {U T M) ■ 5 



(4) 
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By Lemma C.6 and the fact \\U T M(ei - e})|| 2 > V2a k (U T M), we have Pr[E 2 \Ei] > 1 - <5, and thus 
Prf-Ei n E2] > (1 — 25) (1 — 8) > 1 — 35. So henceforth condition on this joint event E\ (~) E 2 . 

Let £r> := II airs ~ _airs|| 2 ._ eo_ anc [ £ ._ — ^£0 ^_ -p^g conditions on iV and the 

bound in (2) implies that £0 < 1 m < 5, so Lemma C.l implies that ah(U T M) > ^/l — e\-a k {M), 
k(U t M) < . l|A f 112 , and that {/ T PairsV is invertible. By Lemma 2.2, 

B[rf) := (L> T Triples(r/)F)(L> T PairsV r )" 1 = (U T M) diag(M T ff)(U T M)' 1 . 
Thus, Lemma C.2 implies 

N A/^ 6/-M1 ^ 1 1 Triples (r/) - Triples (r/)|| 2 | £2 ,_, 

"^ " mh ~ d-eo) -.(Pairs) + ^Pa^s)" (5) 

Let i? := ^Mdiagdl^MeillaJI^M^Ib^.-JI^MefcHs)- 1 and £3 := WEMhMslk^M . Note 
that i? has unit norm columns, and that R~ 1 B(fJ)R = diag(M T fj) . By Lemma C.5 and the fact 
||M|| 2 <>/fc||Af||i = Vfc, 

\\r-% < K (u-M) < —JMh < f (6) 

and 

«W < n^Mf < {1 _ £i) k ak(M yr (7) 

The conditions on TV and the bounds in (2), (3), (4), (5), and (7) imply that £3 < \. By Lemma C.3, 
there exists a permutation r on [k] such that, for all j G [fe], 

\\ s j^j ~ U T ft T (j)/cjh = hjij - R ^T(j)h < 4fc • Il# -1 ll2 • £3 (8) 

where Sj := sign((£j, U T fl T {j))) an d c 'j := II^ T /V(j)ll2 < ||a*t(j)II 2 (*^ e eigenvectors £j are unique up 
to sign Sj because each eigenvalue has geometric multiplicity 1). Since jl T <j) €. range(C/), Lemma C.l 
and the bounds in (8) and (6) imply 

\\SjU£j - Pr^/Cjh < \]\\Sj£j ~ tf T fX r{j )/d j \\l + Wflr^/c'jWl ■ e\ 

< \\sjij - U T f2 T ( j - ) /c , j \\ 2 + \\U T (j) / Cjh ■ £ i 

< 4k- \\R~ l h -e 3 + ei 

< 4fc • . e 3 + £i- 



Therefore, for Cj := Sjc'A1,U£,j), we have 



\ Cj flj - /X T0 -)|| 2 = Wc'jSjUij - p r{j) \\ 2 < \\Pr(j)h ■ (4fc • — - • £3 + El I • 



v/T^ef ■ cr fc (M) 
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Making all of the substitutions into the above bound gives 



\CjHj- H T {j)\\2 



< 



4A: 



1.5 



k 



ek 



flrU)h ~ y/T^e* ■ <r k (M) (1 - el ) ■ a k {MY ^2(1 - e{) ■ a k (M) ■ 5 



+ 



< C- 



/||Triples(r/) — Triples (r/)|| 2 
V (! - e o) • cr fc (Pairs) 

1 1 Pairs — Pairs||2 

l-£o) -o-fe (Pairs) 
k 5 

a k {MY ■ a k (Pairs) 2 • 5 



+ 



| Pairs — Pairs || 2 



1 - e\) • (1 - e - el) ■ <r fc (Pairs) 1 



ln(l/<5) 



N 



A. 3 Some illustrative empirical results 



As a demonstration of feasibility, we applied a modified version of Algorithm A to a subset of arti- 
cles from the "20 Newsgroups" dataset, specifically those in comp . graphics, rec . sport .baseball, 
sci. crypt, and soc. religion, christian, where x\,X2-,%z represent three words from the begin- 
ning (first third), middle (middle third), and end (last third) of an article. We used k = 25 
(although results were similar for k 6 {10, 15, 20, 25, 30}) and d = 5441 (after removing a standard 
set of 524 stop-words and applying Porter stemming) . Instead of using a single if and extracting 
all eigenvectors of B(ff), we extracted a single eigenvector £ x from B{e x ) for several words x £ [d] 
(these x's were chosen using an automatic heuristic based on their statistical leverage scores in 
Pairs). Below, for each such (B(e x ),^ x ), we report the top 15 words y ordered by e^US,x value. 



^V^formatJ 


-B(e'god) 


B(e kcy ) 


-f-H ^polygon/ 


-£> (.CtcamJ 


-B(etoday) 


source 


god 


key 


polygon 


win 


game 


find 


write 


bit 


time 


game 


tiger 


post 


jesus 


chip 


save 


run 


bit 


image 


christian 


system 


refer 


team 


run 


feal 


christ 


encrypt 


book 


year 


pitch 


intersect 


people 


car 


source 


don 


day 


email 


time 


repository 


man 


watch 


team 


rpi 


apr 


ve 


routine 


good 


true 


time 


sin 


public 


netcom 


score 


lot 


problem 


bible 


escrow 


gif 


yankees 


book 


file 


day 


secure 


record 


pitch 


lost 


program 


church 


make 


subscribe 


start 


Colorado 


gif 


person 


clipper 


change 


bit 


fan 


bit 


book 


write 


algorithm 


time 


apr 


jpeg 


life 


nsa 


scott 


wonder 


watch 



The first and fourth topics appear to be about computer graphics (comp. graphics), the fifth and 
sixth about baseball (rec. sports .baseball), the third about encryption (sci. crypt), and the 
second about Christianity (soc. religion. christian). 

We also remark that Algorithm A can be implemented so that it makes just two passes over 
the training data, and that simple hashing or random projection tricks can reduce the memory 
requirement to 0(k 2 + kd) (i.e., Pairs and Triples never need to be explicitly formed). 
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B Proofs and details from Section 3 

In this section, we provide ommitted proofs and discussion from Section 3. 

B.l Proof of Lemma 3.1 
By conditional independence, 

Pi, 2 = EpE[£i ® x 2 \h}} = E\E[xi\h]®E[£ 2 \h]] 

= E[(Mi4) <g> (M 2 e h )} = MAj^wtet ® e t ) M 2 T = M x diag(«;)M 2 T 



t=i 



Similarly, 



A,2,3(rf) = E[E[(xi ® x 2 )(rf,x 3 )|/i]] = E[E[2i|/i] ® E[£ 2 |/i]<»f,E[£ 3 | 

= E[(Mi4) (M 2 e h )(ff, M 3 e h )] = M x [^w t e h ® e^rf, M 3 4) ) M 2 7 



t=i 



Mi diag(M 3 T r/) diag(tt;)M 2 T . 



v-l 



B.2 Proof of Lemma 3.2 

We have U\P\, 2 U 2 = (U^[M\) diag(i?)(M 2 r C/ 2 ) by Lemma 3.1, which is invertible by the assumptions 
on U v and Condition 3.1. Moreover, also by Lemma 3.1, 

Bi,2, 3 (ff) = (U?P lt2:3 (ff)U 2 ) (U^P ly2 U 2 )- 1 

= {U^M 1 diag{M 3 T ff)die,g(w)M^U 2 ) (U^P^U^ 1 

= (^ 1 T M 1 )diag(M 3 T 7?)(^ 1 T M 1 )- 1 {U?M 1 diag(w)M 2 T U 2 ) (U{ P h2 U 2 y l 

= ([/ 1 T M 1 )diag(M 3 T 7?)(^ 1 T M 1 )- 1 {U^P lj2 U 2 ) {U^P^U, 

= (U?M 1 )dmg(M 3 T fj)(U?M 1 )- 1 . 

B.3 Proof of Lemma 3.3 

By Lemma 3.2, 

(U^M 1 )- 1 B 1:2t3 (U 3 e i )(U^M 1 ) = diag(M 3 T L^) 

= diag((4 U 3 T M 3 ^), (9i, *7 3 T M 3 e 2 >, 
= diag(Aj ) i,Ai i 2,. --,Xi,k) 

for all i £ [k], and therefore 



{9 h U 3 T M 3 S k )) 



L 



'(OuUjMae!) (9 ll U 3 T M 3 e 2 ) ••• {9 Xl U 3 T M 3 e k )' 
(4^ 3 T M 3 ei) (9 2 ,U 3 T M 3 e 2 ) ••• (9 2 ,UJM 3 e 3 ) 

{9 k ,U^M 3 e{) (9 k ,U 3 T M 3 e 2 ) ••• (9 k ,U 3 r M 3 e k ). 



6[/ 3 T M 3 . 
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B.4 Ordering issues 

Although Algorithm B only explicitly yields estimates for M3, it can easily be applied to estimate 
M v for all other views v. The main caveat is that the estimators may not yield the same ordering 
of the columns, due to the unspecified order of the eigenvectors obtained in the third step of the 
method, and therefore some care is needed to obtain a consistent ordering. However, this ordering 
issue can be handled by exploiting consistency across the multiple views. 

The first step is to perform the estimation of M3 using Algorithm B as is. Then, to estimate 
Mi-, one mav re-use the eigenvectors in R\ to diagonalize #1,3,2(77% as #1,2,3 (77) and #1,3,2(77) share 
the same eigenvectors. The same goes for estimating M v for other all other views v except v = 1. 

It remains to provide a way to estimate Mi . Observe that M2 can be estimated in at least two 
ways: via the operators #1,3,2(7?), or y ia the operators #3,1,2(77). This is because the eigenvalues 
of #3,1,2(^7) and #1,3,2(7/) are the identical. Because the eigenvalues are also sufficiently separated 
from each other, the eigenvectors #3 of #3,1,2 (if) can be put in the same order as the eigenvectors 
#1 of #1,3,2 by (approximately) matching up their respective corresponding eigenvalues. Finally, 
the appropriately re-ordered eigenvectors #3 can then be used to diagonalize #3,2,1 (jf) to estimate 
Mi. 

B.5 Estimating the mixing weights 

Given the estimate of M3, one can obtain an estimate of w using 

w := M|E[x 3 ] 

where A< denotes the Moore-Penrose pseudoinverse of A (though other generalized inverses may 
work as well), and E[x 3 ] is the empirical average of £3. This estimator is based on the following 
observation: 

E[x 3 ] = E[E[x 3 \h}} = M 3 E[e h ] = M 3 w 

and therefore 

AfjE[£ 3 ] = M\M 3 w = w 

since M3 has full column rank. 

B.6 Proof of Theorem 3.1 

The proof is similar to that of Theorem 2.1, so we just describe the essential differences. As before, 
most perturbation arguments are deferred to Appendix C. 
First, let #1 be the event in which 

||A,2-#l,2||2<Cl,2-/(iM), 
||A,3-#l,3||2<C 1 , 3 -/(^,<5) 

and 

||A,2,3(W " #1,2,3(^)112 < Ci, 2 ,3 • f(N,5/k) 

for all i €. [k]. Therefore by Condition 3.2 and a union bound, we have Pr[#i] > 1 — 35. Second, 
let #2 be the event in which 

,,- *-r , n>, min,-^,-/ \\tlJM 3 (ei — euOlh ■ <5 

7 := minminl^, U^ M 3 (e 3 - e»)| > 3 * 3 " 3 3V J J m 

«6[fcb^/ yek\X)k 
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A max := max \{Bi,U^M^j)\ < jeW " J " (1 + y^H^JS)). 



and 

max je[fc ] \\M 3 ej\ 

i,j£[k] " ' \fk 

Since each 9i is distributed uniformly over S k ~ l , it follows from Lemma C.6 and a union bound 
that Pr[£ 2 |£i] > 1 - 25. Therefore Pr[£i n E 2 ] > (1 - 3<5)(1 — 2(5) > 1 — 55. 

Let C/3 G R dxfc be the matrix of top k orthonormal left singular vectors of M 3 . By Lemma C.l 
and the conditions on N, we have CTk{U 3 U 3 ) > 1/2, and therefore 



min^/ ||M 3 (ei - eV)|| 2 • 5 A max >/efc 3 (l + x /21n(F7^)) ,,», s 

2V^k( k 2 )k 7 " <J 

where 

max,r[ m i II Mae,- II 2 
k'(M 3 ): '-'" " IU 



min^j/ ||M 3 (ej - ej/)|| 2 

Let fji := U^Oi for i G [k]. By Lemma C.l, IJ\P\^2 is invertible, so we may define .61,2,3 (?7i) := 
(^i T i 3 i,2,3(r?i)^2)(^i T Pi,2^2)- 1 . By Lemma 3.2, 

£1,2,3 (#) = (C/ 1 T M 1 )diag(M 3 T ^)(L> 1 T M 1 )- 1 . 

Also define R := UJ M 1 diagdl^Mieil^, \\U± Mx&h, • • • , || ^Afie^ || 2 ) -1 - Using most of the same 
arguments in the proof of Theorem 2.1, we have 

\\R~ 1 \\ 2 <2k(M 1 ), (9) 

k{R) < 4k(Mi) 2 , (10) 



2||A,2, 3 (^) - ^1,2,3(^)112 , 2||Pi j2 , 3 ||2 • ||A,2 " P] 



\Bi, 2 M) - £l,2,3(^)||2 < " ^^- i2±il»£lbf + 



1,2 2 



0-fe(Pl,2) <7fc(P 



1,2 
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By Lemma C.3, the operator Si 2 3(771) has k distinct eigenvalues, and hence its matrix of right 
eigenvectors R± is unique up to column scaling and ordering. This in turn implies that R^ is 
unique up to row scaling and ordering. Therefore, for each i G [k], the Ajj = eJ R\ Bi : 2,3{ffi)Riej 
for j G [k] are uniquely defined up to ordering. Moreover, by Lemma C.4 and the above bounds on 
||-Bi i2)3 (rfj) — B\ ) 2,3(77$) || 2 and 7, there exists a permutation r on [k] such that, for all i,j G [k], 

\ko ~ x i,r(j)\ < (S K (R) + 16k 1 - 5 ■ k(R) ■ WR^Wl ■ A max / 7 ) • \\B lt2 M) - B 1X3 (fJ t ) lb 

< (i2k(M x ) 2 + 256k 1 - 5 ■ k(Mi) 4 • A max / 7 ) ■ ||Bi,2,3(^) - B lx M)h (H) 

where the second inequality uses (9) and (10). Let i)j := (Aij, A 2j -, . . . , Xk,j) G ^ and Vj := 
(Aij, A 2 j, . . . , Afej) G R fc . Observe that Vj = QU 3 M 3 ej = QU 3 T ftz,j by Lemma 3.3. By the 
orthogonality of 0, the fact ||-u|| 2 < vfc||u||oo f° r v G M. , and (11) 

ll 0_1 % - U 3 T p3,TU)h = II ©"Vj - VrV))h 



m - 


-Vr(j)\\2 


Vk- 


W^j ^tU) ii°° 


Vk- 


max Ajj - \i >T (j) 



< (l2y/k • k(M x ) 2 + 256k 2 • «(Mi) 4 • A max / 7 ) ■ WB^M) - ^1,2,3^ 



21 



Finally, by Lemma C.l (as applied to Pi 3 and Pi, 3), 

Up p ll 

II- - II s llr»-l- fl^- 11 , on- 11 Il-TL,3 — -TL,3||2 
llA*3,j - ^3,r(j)\\2 < ||B I/j - U 3 M3,r(i)||2 + 2||// 3)T (j) || 2 , p , • 

Making all of the substitutions into the above bound gives 

11 r, r, \\ < C t s ^V , M/f ^ Mfc/g) f C h2 ,s- f(N,5/k) \\Pi^h-C h2 - f{N/5) 
ll*U " M3,. )l|2 < ? • * • k(Mx) • K (M 3 ) • -^— • ^ -^ + ^—y 2 



C C lt3 -ttN,5) 

6 KJ ' CTk{Pl,3) 



< max II u^ ,-/|| 2 • e. 

j'6[fc] '■"' 

C Perturbation analysis for observable operators 

The following lemma establishes the accuracy of approximating the fundamental subspaces (i.e., 
the row and column spaces) of a matrix X by computing the singular value decomposition of a 
perturbation X of X. 

Lemma C.l. Let X G ]g mxn fr e a matrix of rank k. Let U G ]^j"-xk an( ^ y ^ flj«-xfe fr e matrices 
with orthonormal columns such that range ([/) and range(y) are spanned by, respectively, the left 
and right singular vectors of X corresponding to its k largest singular values. Similarly define 
U G M. mxk and V G R nxk relative to a matrix X G R mxn . Define e x ■= \\X - X\\ 2 , e := ^^, 

and £] : = = £ ° . Assume en < i. Then 

1. ei < 1; 

2. (7 fc (X) = ^((7 T I1>) > (1 - £ ) • <r fc (X) > 0; 



3. a k (U T U)> yr^"'- 



!■ 



I <y k {V T V) > VT^4; 

5. o- k {U^XV)>{l-el)-o- k (X); 

6. for any a G M fc and v G range([/) ; \\Ua — v||| < ||d — ^ T v|| 2 + H^H^ ■ &\- 

Proof. The first claim follows from the assumption on £o- The second claim follows from the assump- 
tions and Weyl's theorem (Lemma E.l). Let the columns of U± G ]R mx ( m - fc ) be an orthonormal 
basis for the orthogonal complement of range([7), so that ||£/j[?7||2 < ex/o~ k {X) < E\ by Wedin's 
theorem (Lemma E.2). The third claim then follows because ||C^ T ^||2 = 1 — ll^i^lll — 1 — £ i- 
The fourth claim is analogous to the third claim, and the fifth claim follows from the third and 
fourth. The sixth claim follows writing v = Ua for some a G R , and using the decomposition 
\\Ua-vf 2 = \\Ua-UU T v\\l+\\U L Ulv\\ 2 2 = \\a-U T v\\l+\\Ul(Ua)\\l < \\a-U T v\\ 2 2 + \\UlU\\ 2 2 \\a\\l < 
|| d— ^7 T w||2 + ||a||2-£i = || d — f^^alli + IMIi'^i where the last inequality follows from the argument 
for the third claim, and the last equality uses the orthonormality of the columns of U. □ 
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The next lemma bounds the error of the observation operator in terms of the errors in estimating 
the second-order and third-order moments. 

Lemma C.2. Consider the setting and definitions from, Lemma C.l, and let Y G B mx " and 

Y G ]R mxn 5e qiven. Define e 2 := 7- — zttt zr and ey := \\Y — Y\\ 2 . Assume En < — ^=. Then 

y J z (l-ej)-(l-eo-ei) * " l|Z u l+v^ 



1. U T XV and U T XV are both invertible, and \\(U r XV)- 1 - (U r XV)- 1 \\ 2 < -^ 



<?k 



TX) : 



2. \\{U-YV){U- XV)- 1 - (U-YV)(U-XV)-% < {l _ £ ^ k{X) + ™f • 



Proof. Let S := U T XV and 6" := J7 T Xy. By Lemma C.l, U T XV is invertible, a k (S) > a k (U T U) ■ 
a k (X) ■ a k {V T V) > (1 - e?) • (T fe (X) (so 5 is also invertible), and \\S - S\\ 2 < £ ■ a k {X) < -^ ■ 

a k (S). The assumption on £0 implies £ % < 1; therefore the Lemma E.4 implies \\SI~ 1 — S" 1 ^ < 
" „ - — l[ 2 /°"fcv -> . — < £2 which proves the first claim. For the second claim, observe that 

l-||5 , -S|| 2 /cr fc (5) cr k (S) - fffc(A)' 

||(L> T yv r )(L> T lv r )- 1 - (c/ T yy)(L> T xy)- 1 || 2 

< ||(t/' T yv r )(c/' T xv r )" 1 - (L> T yy)(f/ T xy)" 1 ||2 + ||(i7 T yv r )(i7 T .x'V')~ 1 - (f/ T yy)(L> T xy)" 1 || 2 

< \\U T YV - U T YV\\ 2 • ||(t/' T XV')- 1 ||2 + \\U T YV\\ 2 • ||(t/' T XV r )- 1 - (U T XV)- 1 \\ 2 

ey ||y||2 ■ £2 

" (1 - s ) • a k (X) + a fc (X) 

where the first inequality follows from the triangle inequality, the second follows from the sub- 
multiplicative property of the spectral norm, and the last follows from Lemma C.l and the first 
claim. □ 

The following lemma establishes standard eigenvalue and eigenvector perturbation bounds. 

Lemma C.3. Let A G M. kxk be a diagonalizable matrix with k distinct real eigenvalues X±, X 2 , . . . ,X k G 
M. corresponding to the (right) eigenvectors £i,£ 2 , ■ ■ ■ ,£& G R all normalized to have \\£i\\ 2 = 1. Let 
R G M. kxk be the matrix whose i-th column is £j. Let A G M. kxk be a matrix. Define €a '■= \\A — A\\ 2 , 
7A := miii^ |Aj — Xj\, and £3 := K( ' €A . Assume £3 < |. Then there exists a permutation r on 
[k] such that the following holds: 

1. A has k distinct real eigenvalues Ai, X 2 , . . . , X k G R, and \X T a\ — Aj| < £3 • 7^4 for all i G [k]; 

2. A has corresponding (right) eigenvectors £1, £2, ••-,£& S R fc , normalized to have \\£,i\\ 2 = 1, 
which satisfy ||£ r (j) — &H2 < 4(fc — 1) • ||-R _1 ||2 ■ £3 /or a// i G [fc]; 

5. £/ie matrix R G R fcxfc whose i-th column is ^ T u\ satisfies \\R — R\\ 2 < \\R — R\\p < Ak 1 ' 2 {k — 
1) ■ \\R-% ■ £3. 

Proof The Bauer-Fike theorem (Lemma E.3) implies that for every eigenvalue Xi of A, there exists 
an eigenvalue Aj of A such that \X\ — Xj\ < ||i? _1 (A — -A)i?||2 < £3 • 7a- Therefore, the assumption 
on £3 implies that there exists a permutation r such that \X T u\ — Aj| < £3 • 7a < -^"- In particular, 



x t - 2 ,x t + 2 



n{Ai,A 2 ,...,A fe } -1, ViG [k]. (12) 
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Since A is real, all non-real eigenvalues of A must come in conjugate pairs; so the existence of a 
non-real eigenvalue of A would contradict (12). This proves the first claim. 

For the second claim, assume for notational simplicity that the permutation r is the identity 
permutation. Let R 6 R fcxfc be the matrix whose i-th column is £j. Define (J & R fe to be the 
i-th row of Rr l (i.e., the i-th left eigenvector of A), and similarly define Q G R fc to be the i-th 
row of RT 1 . Fix a particular i S [k]. Since {£1,^2, ■ ■ ■ ,£k} forms a basis for R , we can write 
£« = ^2j=i c i,j^j f° r some coefficients (H,i,(H,2,--- ,(H,k £ ^- We may assume c^j > (or else we 
replace £j with — &). The fact that ||£j||2 = HCjIh = 1 for all j 6 [k] and the triangle inequality 
imply 1 = || ^ || 2 < Cj^H&l^ + Ej^i \°i,j I II & II 2 = (H,i + Y,&i Kjl> and therefore 

|||i -&H2 < |1 -Cj,j|||^||2 + ^1^11^112 < 2^|qj| 

again by the triangle inequality. Therefore, it suffices to show \cij\ < 2||i? _1 H2 ■ £3 for j 7^ i to prove 
the second claim. 

Observe that A£j = ^4Q^j'=i c i,i'd') = Xw'=i c i,i' A i'£,i'-, and therefore 

fe fc 

^ Ci,i'\>Z,i' + (A- A)ii = A£i = Xiii = Xi ^2 c i,i'd> + (A; - Ai)&. 

Multiplying through the above equation by Q , and using the fact that C/£i' = Mi = *'} gives 

c i,j x j + CI (A ~ A)£i = Xic id + {Xi - Xi)(j£i. 
The above equation rearranges to (Xj — Aj)cjj = (Xi — \i)CJ€i + CJ(A — A)^ and therefore 

1 . IIOIl2-(|A,-A t | + ||(i-A)| t || 2 ) lli?- 1 ^ • (|A, - Xi\ + \\A - A\\ 2 ) 



\Aj — Ai\ \Aj — Ai\ 

by the Cauchy-Schwarz and triangle inequalities and the sub-multiplicative property of the spectral 
norm. The bound \cij\ < 2||i?~ 1 ||2 • £3 then follows from the first claim. 

The third claim follows from standard comparisons of matrix norms. □ 

The next lemma gives perturbation bounds for estimating the eigenvalues of simultaneously di- 
agonalizable matrices Ai, A2, . . . , A^. The eigenvectors R are taken from a perturbation of the first 
matrix A\, and are then subsequently used to approximately diagonalize the perturbations of the 
remaining matrices A2, . . . , A^. In practice, one may use Jacobi-like procedures to approximately 
solve the joint eigenvalue problem. 

Lemma C.4. Let A\,A 2 , ■ ■ ■ ,Ak £ R be diagonalizable matrices that are diagonalized by the 
same matrix invertible R € R fcxfc with unit length columns \\RejW2 = 1, such that each A4 has k 
distinct real eigenvalues: 

R' 1 A i R = diag(Ai,i, Ai )2 , • • • , A ijfc ). 



Let A±,A 2 , . . . , Ak G M. be given. Define ca '■= maxj \\Ai — ^4j||2> 7A := minj min^/ |Ajj — Ajj 
Amax := maxj ,- |Aj ,-|, £3 := . , and £4 := 4/c L5 • ||i? _1 ||2 • £3. Assume £3 < i and £4 < 1. Th 
there exists a permutation r on [k] such that the following holds. 
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1. The matrix A\ has k distinct real eigenvalues \\ t i, Ai^, ■ ■ ■ , Ai^ G R, and \X\j — Ai T (,-)| < 

£3 • 1A for all j G [k] . 

2. There exists a matrix R G M. kxk whose j-th column is a right eigenvector corresponding to 
Aij, scaled so \\RejW2 = 1 for all j G [k], such that \\R — -R r ||2 < WR-iw > w here R T is the 
matrix obtained by permuting the columns of R with r. 

3. The matrix R is invertible and its inverse satisfies ||i? _1 — -R^ 1 !^ < ||-R -1 ||2 ■ £4 • 



1-E4 : 



4- For alii G {2,3, ... ,k} and all j G [k], the (j,j)-th element of R AiR, denoted by A 
e] RT 1 AiRej, satisfies 

£4 \ / £4 



i-:i 



,t>\ f * 1 1 £4 \ 

+ K{R) \— A + 7k^(R) + Tk-—J' £A ' w 

// £4 < 3, i/ten \Xij - \i, T (j)\ < 3£3 • 7A + 4re(.R) • £4 • A max . 

Proof. The first and second claims follow from applying Lemma C.3 to A\ and A\. The third claim 
follows from applying Lemma E.4 to R and R T . To prove the last claim, first define Q G M fc (C7) to 

be the j-th row of R^ 1 (R^ 1 ), and G M fc (0) to be the j-th column of R T (R), so (JA^j = \ >T (j) 
and (J Ai£j = e^ RT 1 AiRej = X\j. By the triangle and Cauchy-Schwarz inequalities and the 
sub-multiplicative property of the spectral norm, 

= \<jM J -tjA i £ j \ 

= \C]{A t - Ai)ij + C]{A % - AiXij - g) + (0 - Cj) T {Ai - AAgj 

+ (0 - 6) T (i* - A i)(ij - S) + (0 - CjVMj + CjMij - S) + (0 - 6) T A(6 - 6)1 

< K7(A - Ai)ij\ + \c]{A t - Ai)(ij - 6)1 + 1(0 - 0) T (A - ^)6l 

+ ko - o) T (A - ^)(6 - 6)1 + ko - 6) t aoi + io t a(6 - 6)1 + ko - oo T ^(6 - 6)1 

< HOH2 • ||A - i^Ha • HOII2 + IIOII2 • ||ii - i^Ha • ||6 - 01b + 110 " 6II2 • \\Ai - AiWASjh 

+ iiCj-6i|2-ii^»-ai|2-ii6-Si|2 

+ 110 - 01b • \\Xi t r(j)Cjh + l|Aj jT (j)0ll2 • 110 ~€jh + 110 -Olb • IIAII2 • 110 -01b- ( 13 ) 

Observe that ||g|| 2 < \\R- l \\ 2 , HOlb < \\Rh, 110 " Olb < \\R~ l - K l h < WR^h ■ j^, 
110 ~~ OII2 ^ 4fc • ll-R -1 ]^ • £3 (by Lemma C.3), and HAII2 < ll-^lb ■ (maxj |Ajj|) • ||i? _1 ||2. Therefore, 
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continuing from (13), |Ajj — \ jT (j)\ ^ s bounded as 



111 £ 4 



\hj ~ \,r(j)\ ^ W R % • \\ R h ■ £A + WR^h • £A • 4/c • \\R- l \\ 2 ■ £3 + ||i2 _1 ||2 ■ ~, e A ■ \\R 

1 — 84. 



+ HiT 1 l 



£4 



2 ' 



1 — £4 



•6A-4/C- p 1 \\ 2 -e 3 



-1- X II f?-l|| £4 

"T A max • it 2 " 



1 — £4 



l-RI 2 + A max ■ lli? IU ■ 4:k ■ \\R H2 • £3 



+ WR- 1 



£4 



2 - 



R\U ■ A max • ll-R lb -4fc-||i? ||2 • £3 



1 — £4 
£3 • 7A + -7= 7— - • £3 • 1A + z £3 • 1A 



y/k-n(R) 



1 — £4 



+ 



£4 



£4 



Vk • k(R) 1 - ^4 



• £3 • 1A 



,„, 1 1 . k(-R) £4 

+ k(K) ■ - • £4 • A max H ^= • £4 • A max H -=- • - • £4 • A max . 

J- — £4 Vk Vk 1 — £4 



Rearranging gives the claimed inequality. 



□ 



Lemma C.5. Let V G R be an invertible matrix, and let R E Mr be the matrix whose j-th 
column is Vej/||V^|| 2 . ZTien ||U|| 2 < k(V), ||-R _1 || 2 < k(V), and «(#) < k(1/) 2 . 

Proof. We have i? = Vdiag(||Vei||2, || Ve2 1|2 3 ■ • • j ||^e*A;||2) _1 , so by the sub- multiplicative property 
of the spectral norm, ||i?||2 < H^lh/minj ||V^-||2 < || V"||2/cfc(V') = n(V). Similarly, ||i? _1 ||2 < 
\\V- 1 \\2-m&x j \\Ve j \\2<\\V- 1 \\2-\\V\\2 = K(V). D 

The next lemma shows that randomly projecting a collection of vectors to M does not collapse 
any two too close together, nor does it send any of them too far away from zero. 

Lemma C.6. Fix any 5 G (0, 1) and matrix A G M mxn (with m < n). Let 9 G W 11 be a random 
vector distributed uniformly over S" 1 ^ 1 . 

minijLj \\A(ei -e})|| 2 • <T 



1. Pr 



mm\{0,A(ei-ej))\ > 



'em 



> 1-5. 



2. Pr 



viG h, \(e,Aii)\< 



\Ae 



i\\2 



III 



1 + y/2\D.{m/5) 



>\-S. 



Proof. For the first claim, let 8q := $/(?) ■ By Lemma F.2, for any fixed pair {i,j} C [n] and 
/3 := Sq/V^, 



Pr 



\(9,A(e t - ej)}\ < \\A(e % - ^)\\ 2 ■-;= — 

■urn \/e 



<exp(^(l-(^/e)+ln(«5o/e))) < <h- 



Therefore the first claim follows by a union bound over all (Tj pairs {i,j}- 

For the second claim, apply Lemma F.2 with /3 := 1 + 1 and t := y21n(m/<5) to obtain 



Pr 



_ 1 1 A p. 1 1 „ 

\(9,Ae i )\> l -^M-(l + t) 



in 



<exp( ^(l-(l + i) 2 + 21n(l+i) 



< exp ( ^ ( 1 - (1 + t) 2 + 2t 



-t 2 /2 



5/m. 
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Therefore the second claim follows by taking a union bound over all i G [m]. □ 

D Proofs and details from Section 4 

In this section, we provide ommitted proofs and details from Section 4. 

D.l Proof of Proposition 4.1 

As in the proof of Lemma 3.1, it is easy to show that 

Qi, 2 ,3($,$) =E[E[x 1 \h]®E[x 2 \h]{$,E[x 3 ® x 3 \h}$}} 

= MiE[4 ® 4(0, (ji 3th <g> fa jh + U 3 ^)]M^ 
= Mi diag«<£ M3,t)(^, /fc.t) + (<?, £s,t$ : * G [fc]) diag(uJ)M 2 T . 
The claim then follows from the same arguments used in the proof of Lemma 3.2. 

D.2 Proof of Proposition 4.2 

The conditional independence properties follow from the HMM conditional independence assump- 
tions. To check the parameters, observe first that 

DI . -I. -1 PT[h2=j\hl=i]-Pl[hl=i] T jyi TTi _ f-.T.rpTj. frp^-l- 

Pr[/n = i\h 2 = j] = p , h - ., = jj^r- = ei diag(vr)r diag(T7r) e,- 

by Bayes' rule. Therefore 

Mie} = E[xi\h2 = j] = OE[e hl \h 2 = j] = Odiag(7r)T T diag(T7?) _1 ej. 
The rest of the parameters are similar to verify. 

D.3 Learning mixtures of product distributions 

In this section, we show how to use Algorithm B with mixtures of product distributions in R n that 
satisfy an incoherence condition on the means fii, fa, ■ ■ ■ , jlk £ ^ n of k component distributions. 
Note that product distributions are just a special case of the more general class of multi-view 
distributions, which are directly handled by Algorithm B. 

The basic idea is to randomly partition the coordinates into £ > 3 "views" , each of roughly the 
same dimension. Under the assumption that the component distributions are product distributions, 
the multi-view assumption is satisfied. What remains to be checked is that the non-degeneracy 
condition (Condition 3.1) is satisfied. Theorem D.l (below) shows that it suffices that the original 
matrix of component means have rank k and satisfy the following incoherence condition. 

Condition D.l (Incoherence condition). Let 5 G (0, 1), £ G [n], and M = [fa\fa\ • • • \fa] G R nxfc 
be given; let M = USV T be the thin singular value decomposition of M, where U G IR ra is a 
matrix of orthonormal columns, S = diag(<7i(M),<72(M), . . . ,<7fc(M)) G R fcxfc , and V G R kxk is 
orthogonal; and let 



n 



c M ■= max - • \\U e} 2 
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The following inequality holds: 

9 n 

cm < ™ 



32 kllnf 

Note that cm is always in the interval [l,n//c]; it is smallest when the left singular vectors in 
U have ±l/y/n entries (as in a Hadamard basis), and largest when the singular vectors are the 
coordinate axes. Roughly speaking, the incoherence condition requires that the non-degeneracy of 
a matrix M be witnessed by many vertical blocks of M. When the condition is satisfied, then with 
high probability, a random partitioning of the coordinates into £ groups induces a block partitioning 
of M into £ matrices Mi, M 2 , . . . , Mi (with roughly equal number of rows) such that the fc-th largest 
singular value of M v is not much smaller than that of M (for each v £ [£}). 

Chaudhuri and Rao (2008) show that under a similar condition (which they call a spreading 
condition), a random partitioning of the coordinates into two "views" preserves the separation 
between the means of k component distributions. They then follow this preprocessing with a 
projection based on the correlations across the two views (similar to CCA). However, their overall 
algorithm requires a minimum separation condition on the means of the component distributions. 
In contrast, Algorithm B does not require a minimum separation condition at all in this setting. 

Theorem D.l. Assume Condition D.l holds. Independently put each coordinate i G [n] into one 
°f £ different sets Ii,I<2, ...,!(. chosen uniformly at random. With probability at least 1 — 5, for 
each v G [£], the matrix M v G M) Xv ' xk formed by selecting the rows of M indexed by X v , satisfies 

(?k{M v ) > a k (M)/(2V~£). 

Proof. Follows from Lemma D.l (below) together with a union bound. □ 

Lemma D.l. Assume Condition D.l holds. Consider a random submatrix M of M obtained by 
independently deciding to include each row of M with probability l/£. Then 



Pr 



o- k {M) > a fe (M)/(2v^)l > 1 - 5/£. 



Proof. Let z±,Z2, ■ ■ ■ ,z n £ {0,1} be independent indicator random variables, each with Pr[zj = 
1] = l/£. Note that M T M = M T diag(zi, z 2 , • • • , z n )M = Ya=i ZiM T aejM, and that 

a k (Mf = A min (M T M) > X min (S) 2 ■ X m Jj2z t U T e t eJu\ 

Moreover, < z^aejU < (k/n)c M I and A min (EE" =1 Zi U T eiejU]) = l/£. By Lemma F.3 (a 
Chernoff bound on extremal eigenvalues of random symmetric matrices) , 



Pr 



V/=i / 



< k . e -(3/4)V(2fc M fc/n) < S/£ 



by the assumption on cm- D 
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D.4 Empirical moments for multi-view mixtures of subgaussian distributions 

The required concentration behavior of the empirical moments used by Algorithm B can be easily 
established for multi-view Gaussian mixture models using known techniques (Chaudhuri et al., 
2009). This is clear for the second-order statistics P a j, for {a,b} G {{1, 2}, {1, 3}}, and remains 
true for the third-order statistics A,2,3 because X3 is conditionally independent of x\ and x 2 given 
h. The magnitude of {U^9i,x^) can be bounded for all samples (with a union bound; recall that 
we make the simplifying assumption that P\^ is independent of -Pi, 2, 3, and therefore so are U3 and 
-Pi, 2, 3)- Therefore, one effectively only needs spectral norm error bounds for second-order statistics, 
as provided by existing techniques. 

Indeed, it is possible to establish Condition 3.2 in the case where the conditional distribution 
of x v given h (for each view v ) is subgaussian. Specifically, we assume that there exists some a > 
such that for each view v and each component j G [k], 



E 



expf X(u,cav(x v \h = j) 1/2 (x v -E[x v \h = j]))j 



<exp(aA 2 /2), VAgM,uGS 



d-l 



where cov(x|/i = j) := E[(x„ — E[x^|/i = j]) <g> (x v — E,[x v \h = j])\h = j] is assumed to be positive 
definite. Using standard techniques (e.g., Vershynin (2012)), Condition 3.2 can be shown to hold 
under the above conditions with the following parameters (for some universal constant c > 0): 



H'n 



mm w-j 

je[k] 



AT orV^ + logCL/^ a 3 / 2 (d + log(l/<5)) 
JVQ := log 



W'r, 



Wr, 



C a)b := c- (maxj||cov(x„|/i = j)\\l /2 , ||E[x„|/t = j]|| 2 : v G {a,b},j G [k]\j 
Ci,2,3 := c- fmax{||cov(x„|/i = j)\\l' 2 , \\E[x v \h = j}\\ 2 ■ v G [3],j G [k]\) 



/•( V. S) : J kn °^ l l 5 ) 1 > l<* 3 ' 2 y/log(N/5)(d + log(l/<5)) 



N 



W m \nN 



E General results from matrix perturbation theory 

The lemmas in this section are standard results from matrix perturbation theory, taken from Stewart 
and Sun (1990). 



Lemma E.l (Weyl's theorem). Let A, E G 



with m > n be given. Then 



max\ai(A + E) - <Ji(A)\ < \\E\\ 2 . 

ie [n] 

Proof. See Theorem 4.11, p. 204 in Stewart and Sun (1990). □ 

i mxn with m> n be given. Let A have the singular 



Lemma E.2 (Wedin's theorem). Let A,E G 
value decomposition 



u 2 



3 J 



A[V l V 2 



Si 

s 2 
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Let A := A + E, with analogous singular value decomposition (U±, U2, U3, Ei,S2, V1V2). Let <& be 
the matrix of canonical angles between range(£/i) and range (C/i), and Q be the matrix of canonical 
angles between range(Vi) and range(Vi). // there exists S, a > such that minjcrj(Si) > a + 5 and 
maxj<Ti(X2) < a, then 

Pj 9 

max{|| sin$||2, || sinSj^} < ~ 



5 



Proof. See Theorem 4.4, p. 262 in Stewart and Sun (1990). 



□ 



Lemma E.3 (Bauer-Fike theorem). Let A, E £ R kxk be given. If A = V diag(Ai, A2, • • • , \k)V 1 
for some invertible V £ M. , and A := A + E has eigenvalues Ai, A2, • • • , A&, then 



maxmin|Aj — AJ < II V EV\ 
ie[k]je[k]' 

Proof. See Theorem 3.3, p. 192 in Stewart and Sun (1990). 



2- 



D 



Lemma E.4. Let A,E G M. kxk be given. If A is invertible, and \\A l E\\2 < 1, then A := A + E is 

invertible, and 

H Pll II i-ill 2 

1 -C/ o /i o 



li- 1 -^" 



< 



Proof. See Theorem 2.5, p. 118 in Stewart and Sun (1990). 



□ 



F Probability inequalities 

Lemma F.l (Accuracy of empirical probabilities). Fix ft = (/j,i,/j,2, ■ ■ ■ ,Hn) G A m_1 . Let x be a 
random vector for which Pv[x = e\\ = ^ for all i G [m], and let x\,X2, ■ ■ ■ ,x n be n independent 
copies of x. Set fi := (1/n) Yl7=i %i- F° r °^ t > 0, 



Pr 



I A* — A*|| 2 > 



l + y/i 



n 



< e 



Proof. This is a standard application of McDiarmid's inequality (using the fact that ||/x — fx\2 has 
\J2/n bounded differences when a single Xi is changed), together with the bound E[||/i — /X| 1 2] < 
1/y/n. See Proposition 19 in Hsu et al. (2012). □ 

Lemma F.2 (Random projection). Let 6 G W 1 be a random vector distributed uniformly over 
iS" -1 , and fix a vector v G R n . 



1. If f3 G (0,1), then 



Pr 



1(0,0)1 <N| 2 - — •/? 

n 



<exp(^(l-/3 2 + ln/3 2 ; 



2. If/3> 1, then 



Pr 



|(0,0)| > ||0|| 2 - — ■/? 



» 



<exp(^(l-/3 2 + ln/3 2 )). 



Proof. This is a special case of Lemma 2.2 from Dasgupta and Gupta (2003). 



□ 
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Lemma F.3 (Matrix Chernoff bound). Let X±, X2, ■ ■ ■ , X n be independent and symmetric m x m 
random matrices such that ■< Xi -< rl, and set I := A m i n (E[Xi + X2 + • • • + -X'n])- F° r an V 
£€[0,1], 



Pr 



A, 






e)-l 



<m . e -r»I/(2r)_ 



Proof. This is a direct corollary of Theorem 19 from Ahlswede and Winter (2002). 



□ 



G Insufficiency of second-order moments 

Chang (1996) shows that a simple class of Markov models used in mathematical phylogenetics 
cannot be identified from pair-wise probabilities alone. Below, we restate (a specialization of) this 
result in terms of the document topic model from Section 2.1. 

Proposition G.l (Chang, 1996). Consider the model from Section 2.1 on (h,Xi,X2, ■ ■ ■ ,xe) with 
parameters M and w. Let Q G R be an invertible matrix such that the following hold: 

1. l T Q = l T ; 

2. MQ^ 1 , Q diag(u>)M T diag(MtZfp 1 , and Qw have non-negative entries; 

3. Q di&g(w)Q T is a diagonal matrix. 

Then the marginal distribution over (xi,X2) is identical to that in the case where the model has 
parameters M := MQ~ 1 and w := Qw. 

A simple example for d = k = 2 can be obtained from 



M :-- 



p 1 — p 
1 — p p 



w :- 



1/2' 

1/2 



Q:= 



l+ v / l+4p(l-p) 
p 2 

_ i-yi+4 P (i-p) 

j. y 2 



for some p G (0, 1). We take p = 0.25, in which case Q satisfies the conditions of Proposition G.l, 
and 



M 

M = MQ- 1 



0.25 0.75 
0.75 0.25 

0.6614 0.1129 
0.3386 0.8871 



w 



w 



)w 



"0.5" 




0.5 


•) 


"0.7057" 


0.21 


)43 



0.3125 0.1875 
0.1875 0.3125 



0.2188 0.0938 
0.0938 0.0938 



In this case, both (M, w) and (M, w) give rise to the same pair- wise probabilities 

Mdiag(w)M T = Mdiag^M 1 
However, the triple-wise probabilities, for r] = (1,0), differ: for (M,w), we have 

M diag(M T ??) diag(w)M T 

while for (M, w), we have 

ma- ni>T \a- t-MT/rT [0-2046 0.1079 
Mdiag(M r 1 )dmg{w)M « QmQ QmQ6 
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